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ABSTRACT 

We investigate the clustering properties of ^1550 broad-line active galactic nuclei (AGNs) at (z)=0.25 
detected in the ROSAT All-Sky Survey (RASS) through their measured cross-correlation function 
with ^46000 Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey. By measuring the 
cross-correlation of our AGN sample with a larger tracer set of LRGs, we both minimize shot noise 
errors due to the relatively small AGN sample size and avoid systematic errors due to the spatially- 
varying Galactic absorption that would affect direct measurements of the auto-correlation function 
(ACF) of the AGN sample. The measured ACF correlation length for the total RASS- AGN sample 
(< L^- 2A ^ V >= 1.5 x 10 44 ergs" 1 ) is r = 4.3±°; 4 h~ x Mpc and the slope 7 = 1.7±°' 4 . Splitting the 
sample into low and high Lx samples at £!^ 5 ~ 10kcV = 10 44 ergs -1 , we detect an X-ray luminosity- 
dependence of the clustering amplitude at the ~2.5<r level. The low Lx sample has ro = 3.3lg g 
h~ x Mpc (7 = 1.71q'|), which is similar to the correlation length of blue star-forming galaxies at 

low redshift. The high Lx sample has tq = 5.4j^'o Mpc (7 = 1.91q'2), which is consistent 
with the clustering of red galaxies. From the observed clustering amplitude, we infer that the typ- 
ical dark matter halo (DMH) mass harboring RASS-AGNs with broad optical emission lines is log 

(Mdmh/O -1 M q )) = 12.618 J, H- 8 !^ 6 ) iS-iloi for the total > low L x, and high L x RASS-AGN 
samples, respectively. 

Subject headings: galaxies: active — large-scale structure of universe — X-rays: galaxies 



1. INTRODUCTION 

Galaxies and active galactic nuclei (AGNs) are not dis- 
tributed randomly in the universe. The small primordial 
fluctuations in the matter density field present in the very 
early universe have progressively grown through gravita- 
tional collapse to create the complex network of clusters, 
groups, filaments, and voids seen in the distribution of 
structure today. Galaxies and AGNs, as well as groups 
and clusters of galaxies, are believed to populate the col- 
lapsed dark matter halos (DMHs). The clustering of 
galaxies and AGNs therefore reflects the spatial distribu- 
tion of dark matter in the universe. This allows clustering 
measureme nts to be used to derive cosmological param - 
eters (e.g., iPeacock et al.l [20011 : lAbazajian et al.ll2005l ). 
However, these measurements also allow us to study the 
complex physics which governs the creation and evolu- 
tion of galaxies and AGNs, as well as the co-evolution 
of galaxies and AGNs. The co-evolution scenario is mo- 
tivated by the observed correlation between the mass of 
the central super-massive black hole ( SMBH) and the 
stella r velocity dispersion in the bulge (jGebhardt et al.l 
120001 iFerrarese fc Merritti 120001 ). lending strong evi- 
dence to an interaction or feedback mechanism between 
the SMBH and the host galaxy. The specific form of 
the feedback mechanism, as well as the details of the 
AGN triggering, accretion, and fueling mechanisms, re- 
mains unclear. Different cosmological simulations ad- 
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dress possible scenarios fo r the co-evolution of AGNs and 
their host galaxies (e . g., [K auffmann & Hac hneltl 120001 : 
iDi Matteo et al.l 120051 : ["Cattaneo et al.ll2006f ). Large vol- 
ume, high-resolution simulations that include physical 
prescriptions for galaxy evolution and AGN feedback 
make predictions for the spatial cluster ing and large- 
scale e nvironments of AGNs and ga l axies (ISpringel et all 
[20051: iColberg fc Pi Matted [20081: iBonoli et al l 120091) . 
Observed clustering measurements of AGNs can be used 
to test these theoretical models, put constraints on the 
feedback mechanisms, identify the properties of the AGN 
host galaxies, and understand the accretion processes 
onto SMBHs and their fueling mechanism. 

X-ray surveys allow us to identify AGN activity with- 
out contamination from the emission of the host galaxy, 
i.e., therefore efficiently detecting even low luminosity 
AGNs. In the current era of deep and wide-area X- 
ray surveys with extensive spectroscopic follow-up, mea- 
surements of the three-dimensional (3D) cl ustering of 
AGN in various redshift ranges are emer ging (jCoil et al.l 
120091: iGilli et al.ll2005t lYang et al.ll2006D . However, our 
knowledge of AGN clustering in the low redshift universe 
(z < 0.4) is still poor, except for optically selected type 
II AGNs dWake et al.ll200llLi et all 12009) . This is due 
to the lack of observable comoving volume and the low 
number density of low-z AGNs. Exceptionally large sur- 
vey areas with good sensitivity are needed to acquire a 
sufficiently large number of objects for clustering mea- 
surements. 

To date the ROSAT All-Sky Survey (RASS; 
IVoges et al.l 119991 ) is the most sensitive survey to have 
mapped the entire sky in X-rays. Surveys with mod- 
ern higher-sensitivity X-ray observatories such as XMM- 
Newton and Chandra cover much smaller areas of the 
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sky (area: ~0.1-10 deg 2 ). Therefore, the available co- 
moving volume from these deeper data sets is not suffi- 
cient to accurately measure AGN clustering at low red- 
shifts (z 5, 0-4). Serendipito us surveys, su ch as extended 
ChAMP (jCovev et al.ll200ll and 2XMM (|Watson et alJ 
2009), cover larger areas (~33 deg 2 and ~360 deg 2 , re- 
spectively). However, the large variations in sensitivity 
between different observations and the non-contiguous 
sky coverage make serendipitous surveys unsuitable for 
wide-area clustering measurements. 

A few studies have attempted to measure the auto- 
correlation function (ACF) of RASS-based AGN sam- 
ples. Two-dimensional ( 2D) angular correlation func- 
tions ()Akvlas et al.l l2000f) do not require redshift mea- 
surements for each AGN, but the projection heavily di- 
lutes the clustering signal. Furthermore, the deprojec- 
tion of the angular clustering to the 3D correlation func- 
tion is subject to uncertainties in the redshift distribu- 
tion, which can be substantial. Direct measurements 
of the 3D RASS-AGN ACF with s pectroscopic redshift 
mea surements have been made by iMullis et al.l (|2004T > 
and iGrazian et all ()2004f ) with a few hundred AGNs, re- 
spectively, providing clustering measurements that have 
large statistical uncertainties caused by the relatively 
s mall sample size. 

lAnderson et alJ (|2003l I2007D positionally cross- 
correlated RASS sources with spectroscopic data 
available from the Sloan Digital Sky Survey (SDSS). 
This dramatically increased the number of RASS- AGNs 
with spectroscopic redshift measurements, which we 
use here to provide significant improvements in the 
measurement of AGN clustering at low redshift. Fur- 
thermore, the availability of spectroscopic redshifts for 
large samples of SDSS galaxies in the same volume 
allows us to use an alternative approach to infer the 
clustering of AGN using calculations of the AGN-galaxy 
cross-correlation function (CCF). This approach uses 
much larger samples of AGN-galaxy pairs and hence 
significantly reduces the uncertainties in the spatial 
correlation function compared with direct measurements 
of the AGN ACF. Furthermore, the use of a CCF avoids 
the problem that we have to correct for the complex 
angular dependences of limiting sensitivity in the X-ray 
sample. 

We therefore have initiated a program to investigate 
the clustering properties of low redshift {z ~ 0.25) RASS- 
AGNs through measurements of the CCF of these AGNs 
with SDSS Luminous Red Galaxies (LRGs). In this 
study, we chose LRGs as the corresponding galaxy sam- 
ple because they have a significant overlap in redshift 
range with our X-ray sample (details arc described later). 

In this first paper of a series, we explain the data se- 
lection, as well as the calculation of the CCF and the 
inferred RASS-AGN ACF. We also investigate the X- 
ray luminosity dependence of the clustering properties 
and biases. In a follow-up paper (T. Miyaji et al., in 
preparation), we will focus on applying the halo occu- 
pation distribution (HOD) model to the calculated CCF 
between RASS- AGNs and LRGs. 

The paper is organized as follows. In Section 2, we 
describe the construction and properties of the LRG 
and X-ray AGN samples in details. All essential steps 
to measure the CCF, compute the ACF via the CCF, 
and estimate errors are explained in Section 3. The re- 



sults of the clustering measurements for the different X- 
ray AGN samples and their luminosity dependence are 
given in Section 4. We discuss these results in Sec- 
tion 5 in the context of other studies and conclude in 
Section 6. Throughout the paper, all distances, lumi- 
nosities, and absolute magnitudes are measured in co- 
moving coordinates and given in units of hr x Mpc, where 
h = i/o/100kms _1 , unless otherwise stat ed. We use a 
cosm ology of 17m = 0.3 and fl\ = 0.7 (jSpergel et alJ 

2003) . We use AB magnitudes throughout the paper. 
All uncertainties represent a la (68.3%) confidence in- 
terval unless otherwise stated. 

2. DATA 

2.1. SDSS Luminous Red Galaxy Sample 

The opti cal data analyzed in this study are drawn from 
the SDSS (jYork et al.ll2000D . We use data both from the 
main galaxy sample, which has a spectroscopic depth 
of r = 17.7 ([Strauss et al.ll2002j ). and the LRG sample, 
which has a spectroscopic depth of r = 19.5, significantly 
fainter than the main galaxy sample. The LRG sample 
was designed for studies of large-scale structure to higher 
redshift; it covers a larger volume than the main galaxy 
sample. Here we use the LRG sample as a large-scale 
structure tracer set to calculate the CCF with the RASS- 
AGN, as the LRG sample covers a similar redshift range 
as the RASS-AGNs. 

The target selection and the properties of the LRG 
sample are described in detail in Eiscnstei n et al.l (2001). 
Two different selection criteria ('cut I' and 'cut IF) were 
introduced in identifying LRGs as at z > 0.4 the typical 
4000 A break in the spectral energy distribution (SED) 
of an early- type galaxy drops out of the g band and falls 
into the r band. lEisenstein et alJ (|2001l) shows that up to 
z ~ 0.38 a volume-limited sample of LRGs with passively 
evolving luminosity and rest-frame colors is selected with 
a very high effici e ncy ( 95% for cut I and 90% for cut II). 
lEisenstein et al.l (|2001h advised that the LRG selection 
algorithm should not used for objects z < 0.15. 

2.1.1. Extraction of the SDSS Luminous Red Galaxy 

Sample 

The LRG sample was obtained using the web-based 
SDSS Catalog Archive Server Jobs System 3 . The ap- 
propriate objects are selected through the prime target 
flag 'galaxy _red' (jEisenstein e t al.l l2001ft In SPSS DR2 
the m odel magnitude code was changed (|Abazaiian et al.l 

2004) to improve the star-galaxy separation. This also 
caused a slight change in the LRG sample definition. We 
make use of the updated LRG sample selection. 

We extract 115,577 LRGs in the DR7 data release 
that have a spectroscopic redshift of z > 0.15, a galaxy 
spectral-type classification, and a redshift confidence 
level of Zconf. > 0.95. We discovered that 3% of ob- 
jects flagged 'galaxy _red' (i.e., both cut selections) in 
DR7 do not fall with in the LRG selection criteria (see 
lEisenstein et al.ll2001l ). This does not happen in earlier 
data releases prior to DR7. 

To calculate the RASS-AGN-LRG CCF, we wish to 
define an LRG sample that is both volume-limited and 
contains a high number density of LRGs to maximize 

3 http : //casjobs . sdss . org/CasJobs/ 
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th e number of possib le AGN-LRG pairs. Table 1 
of IZehavi et al.l ()2005al ) defines three different volume- 
limited LRG samples, corresponding to different lumi- 
nosity ranges, that they use to measure the LRG ACF. 
The first volume-limited LRG sample, with —23.2 < 
M g < -21.2 and 0.16 < z < 0.36, contains the high- 
est number density of objects. Here we adopt the same 
LRG sample definition, as this sample is the best suited 
to our scientific goals. We select LRGs from DR7 that 
meet the criteria of this sample, and we refer to the sub- 
sequent data here as 'LRG sample'. We further limited 
our survey area to DR4+ to match the X-ray sample (see 
Section [2~3)) . The properties of this sample are shown in 
Table rjj 

The selected LRGs span a range in redshift. To con- 
struct a volume-limited sample, one must correct their 
time-evolving SED to account for the evolution of the 
stellar populatio n and the redshi f ting. As described in 
the Appendix of lEisenstein et all (|2001h . we use the ex- 
tinction corrected r* etro magnitude to construct the k- 
corrected and passively evolved r est-frame ffp et ro mag- 
nitudes (non-star- forming model). IZehavi et al.1 |2005a) 
assively evolve their M a to z = 0.3 instead of z = 
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Eisenstein et al.ll2001t we refer to this magnitude by us- 
ing the notation Mg =0 - 3 ). The stellar population synthe- 
sis (SPS) code described in IConrov et al.l (|2009l ) is used 
to derive the purely passive evolution correction. LRGs 
are AM g = —0.27 brighter at z = 0.3 comp ared to z = 0. 



This v alue is similar to the one used in IZehavi et al.l 
(2005aj) who state that the evolution is about 1 mag per 
unit redshift. 

We extensively tested that ou r LRG sample m e et the 
LR G selection criteria o f both Eisenste in et al.l (j2001l ) 
and lZehavi et al.l (|2005a| ). We only consider objects that 
are located in areas of the sky that have a spectroscopic 
completeness fraction in DR7 of / comp i > 0.8. The vast 
majority of the objects in our sample (99.9%) belong to 
the LRG 'cut I' criteria. 

2.1.2. Accounting for the SDSS Fiber Collision 

iBlanton et al.l (|2003af) describe the algorithm used for 
positioning tiles in the plane of the sky and assign- 
ing spectroscopic fibers for observing objects in the 
SDSS. The resulting spectroscopic sample of SDSS is 
spatially biased in that fibers on the same tile cannot 
be placed closer than 55 arcsec. Spectroscopic red- 
shifts are available for galaxy pairs closer than 55 arc- 
sec only if the field is observed at least twice. As dis- 
cussed in Zchavi et al. (2002, 2005a), not taking into ac- 
count the effects of these fiber collisions would systemati- 
cally underestim ate measured correlati on functions, even 
on large scales. IBlanton et al.l (|2005af) provide publicly 
available collision-corrected SDSS catalogs suitable for 
robust large-scale structure studies. They assign a red- 
shift to a galaxy which has no spectroscopic redshift by 
giving it the redshift of that galaxy in a galaxy group 
that is positionally the closest and has a spectroscopic 
redshift. However, these catalogs are generated only for 
the main galaxy sample, not the LRG sample. 

To correct our LRG sample for fiber collisions, we 
therefore use the following approach. We select from the 
SDSS catalog Archive Server all LRG objects that passed 
the pure photometric-based cut I and cut II selection 



criteria, with no spectroscopic restrictions applied. Pho- 
tometrically selected LRGs that get a redshift assigned 
are only those which have a spectroscopic LRG within a 
separation of d < 55'.'0. This corresponds to 2% (1004 
objects) of the total LRG sample. Our redshift assign- 
ment is a two stage process. First, if the photometrically 
selected LRG is included in the primary spectroscopic 
LRG sample but its redshift was initially rejected due to 
a low confidence level of z con f. < 0.95, the spectroscopic 
redshift is now accepted independent of its confidence 
level to correct for fiber collisions. This specific proce- 
dure assigns a redshift only to nine objects. 

As a next step, we make use of the substantially im- 
prove d DR7 photometric redshift code (|Abazaiian et al.l 
2009) compared to earlier SDSS data releases. If the dif- 
ference between the spectroscopic redshift (z spcc j) of an 
LRG and the photometric redshift (z p hoto,i) of a neigh- 
boring (within a 55 arcsecond radius) LRG is 



I ^spec,j -£photo,i |^ ^•Zplioto,; 



(1) 



then the photometric LRG is assigned the same redshift 
as the spectroscopic LRG. Otherwise the photometric 
redshift is taken to be the correct redshift for the photo- 
metric LRG. More than half of all photometric LRG red- 
shifts are assigned the spectroscopic redshift of a neigh- 
boring LRG, while the remaining objects are assigned 
their photometric redshift. 

The above procedure for redshift assignment is used 
because simply adopting the photometric redshift for all 
objects that do not have a measured spectroscopic red- 
shift would smear out the clustering signal in redshift 
space and there fore dilute the amplit ude of the corre- 
lation function. IBlanton et al.l (|2003al) verified that for 
areas of the sky with overlapping spectroscopic tiles, 60% 
of the galaxies that had a neighboring galaxy within 
55 arcsec were with in 10 Mpc of the neighbor galaxy. 
IZehavi et al.l (|2005al ) demonstrated that LRGs are highly 
clustered, so using this procedure for LRGs will result in 
a higher success rate. The average photometric redshift 
error [bzj (1 + z)) is 2%, and within our sample less than 
1% (437 objects) of the LRGs are assigned a photometric 
redshift. 

2.2. The X-ray Sample 

The RASS (|Voges et all 119991) is currently the most 
sensitive all-sky survey in the X-ray regime, with a typ- 
ical flux limit of f x ~ 10~ 13 erg cm" 2 s" 1 (0.1-2.4 keV). 
The area of the sky covered is 99.7%, all of which has 
at least 50 s of exposure with ROSAT Position Sensitive 
Proportional Counters (PSPC). 

Based on the SP SS data release 5 (DR5), 
lAnderson et all (|2007l ) classify 6224 AGNs with 
broad permitted emission lines in excess of 1000 kms -1 
FWHM and 515 narrow permitted emission line AGNs 
matching RASS sources within 1 arcmin (note that 
86% of all matches fall within 30'.'0). Consequently, 
broad-line AGNs account for 92% of all RASS/SDSS 
classifications. In this study, we use the broad emission 
line RASS-AGNs only. The RASS sourc es are taken 
from both the RASS Bright source catalog fVo ees et al 



1999) and the RASS Faint source catalog ([Voges et al 

2000) and have a maximum likelihood ML > 7. 4 



1 lAnderson et alj (120031 . 120071 ) say that sources with ML > 10 
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Fig. 1. — 0.1-2.4 keV X-ray luminosity vs. redshift for the 
broad emission line A GN sample in SDSS DR5 determined by 
I Anderson et al.l {2007). Black symbols show objects in the red- 
shift range 0.16 < z < 0.36, which define the AGN sample 
used here. The dashed line corresponds to the common divid- 
ing line between low luminous AGNs and high luminous QSOs 
(LY 5_1 ° keV = 10 44 ergs _1 ) defining the low and high luminosity 
RASS-AGN samples used here. 



The data cover an area of 5740 dee 2 . The derived 
luminosities in Ander son et al.l (|2007 ) are based on a 
flat ACDM cosmology with (fi M , fi A , ff )=(0.3, 0.7, 
70kms- 1 Mpc- 1 ). 

The RASS/SDSS target selection consists of two se- 
lection algorithm: the main SDSS algorithm for galaxies 
and AGNs (15.0 < m < 19.0) and the specific assign- 
ment of S PSS fibers for RASS i dentifications (15.0 < 
m < 20.5). lAnderson et al.l (|2003f ) explain in detail their 
selection of optical counterpart candidates. 

To understand possible bias effects, which may influ- 
ence our CCF measurements, we summarize the selection 
method briefly. The algorithm considers those SDSS op- 
tical objects within 1 arcmin of the X-ray source position. 
Optically detected objects qualify for SDSS spectroscopy 
if their <?, r, or i magnitude is 15.0 < m < 20.5. Since not 
every optically detected object that fulfills the criteria 
can be observed spectroscopically, different priority lev- 
els are introduced in target selection. The highest prior- 
ity is given to objects with a triple positional coincidence 
of RASS X-ray source, SPSS opt ical object, and FIRST 
(|Becker. White, fc HelfandlH99ai radio source. The next 
priority level consists of SPSS objects in the RASS er- 
ror circles with unusual optical colors indicating AGNs 
(e.g., UV-excess with u — g < 0.6). The third priority 
includes less likely, but plausible, counterparts such as 
bright stars and galaxies. The last group consists of any 
object in the magnitude limit that falls within the RASS 
error circle. One RASS/SPSS spectroscopic fiber is as- 
signed to the object with the highest priority level and a 
spatial offset less than 27.5 arcsec from the X-ray source. 
If multiple objects populate the same priority level, one 
is chosen at random for SPSS spectroscopy. 

Before consideration is given to potential RASS tar- 
gets, an SPSS fiber is first assigned to objects that belong 

have been included. However, investigating the sample shows that 
ML > 7 has been used. 



to the main SPSS sample, which includes both galax- 
ies and optically identified AGNs. About 80% of all 
RASS/SPSS identifications are independently targeted 
for spectroscopy by the SPSS main target selection al- 
gorithm. However, objects in the main sample can be 
observed at the cost of excluding high-priority RASS tar- 
gets. In these cases, the RASS/SPSS identification is 
missing. For these reasons, the RASS/SPSS AGN sam- 
ple is only reasonably complete (> 90%) for AGNs with 
15.0 < m < 19.0. 

Objects identified by RASS/SPSS priority levels and 
not because of their belonging to the main SPSS sample 
tend to be optically fainter, on average, than the main 
AGN sample. In the u — g versus g — r color-color space, 
objects from both selection algorithms lie in the same 
region. Studying their redshift distributions shows that 
objects identified by the RASS/SPSS priority levels are, 
on average, at higher redshifts than the identifications re- 
lying on the main SPSS sample. This also explains their 
lower observed count rates X-ray fluxes. The X-ray hard- 
ness ratios of both populations are very similar. Being 
fainter in X-rays and in the optical wavelengths leads to 
the same X-ray to optical flux ratio for the RASS/SPSS 
priority levels AGNs and the main SPSS sample AGNs. 
Therefore, except that RASS/SPSS priority levels AGNs 
are found on average at higher redshifts, no intrinsic dif- 
ferences between the RASS/SPSS AGNs from both se- 
lection methods are detected. 

To create a well-defined X-ray AGN sample to use for 
clustering measurements, we focus solely on broad emis- 
sion line AGNs. ROSATs sensitivity to the soft energy 
band that is limited to <2.4 keV selects against X-ray ab- 
sorbed AGNs (type II), which are usually optically clas- 
sified as narrow lines objects. Therefore, the sample of 
broad emission line RASS-AGN is much more complete 
than the narrow emission line RASS-AGN sample. Fur- 
thermore, the X-ray/optical counterpart identification is 
more reliable for the broad emission line AGNs than for 
the narrow emission line AGNs, which hav e much higher 
surface density (see lAnderson et al.l 2003, 20071 for de- 
tails). 

In addition to studying the clustering of the RASS- 
AGN sample, here we also test for possible differences 
in the clustering of low versus high X-ray luminosity 
RASS-AGN (Figure[J). The commonly used X-ray divid- 
ing line between Seyfert AGNs and high luminous QSOs 
is L^ 5 - wkeV = 10 44 ergs" 1 in the 0.5-1 keV band (in- 
trinsic luminosity; Ma inieri et al.l 120021 ) . Although this 
dividing line is somewhat arbitrary, i t is widely em- 
ployed in literature, and we use it here. lAnderson et al.l 
( 2007) provide a table which list the galactic absorption- 
corrected 0.1-2.4 keV luminosities assuming a photon 
index of T = 2.5. With this index, the AGN/QSO 
dividing line corresponds to a 0.1-2.4 keV luminosity 
of log (Lx/iergs- 1 )) = 44.298 ( Figure ED). This index 
agrees well with that found by iPiconcelli et al.1 (2005) 
in the XMM-Newton spectra of PG quasars in the 0.5- 
2 keV band, where f = 2.73^q In the same pa- 
per, it is shown that at higher energies (2-12 keV) the 
mean photon index for PG quasars gets harder and is 
r = 1.89 ± 0.11. He r e we use T = 2.5 to be consistent 
with lAnderson et al.l (|2007l ). 

The properties of our total RASS-AGN sample, low Lx 
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RASS-AGN sample, and the high L x RASS-AGN sam- 
ple are shown in Table [TJ In order to measure the CCF 
of the RASS X-ray AGNs with LRGs, the X-ray samples 
must cover the same area of sky and redshift range as the 
LRG sample used here. In all samples, no X-ray detected 
AGN is also classified as an LRG. However, we cannot 
exclude the possibility that RASS-AGNs are hosted by 
LRGs and outshine their host galaxies. In contrast to the 
high L x RASS-AGN sample, the total RASS-AGN sam- 
ple and the low Lx RASS-AGN sample are not volume- 
limited. We calculate the comoving number density in 
Table [1] of these two samples in the following way. For 
a specific R.A. and decl. (contained in the DR4+ ge- 
ometry), we determine the galactic absorption value Nn 
and the RASS exposure time. The RASS faint source 
catalog contains sources with at least six source counts. 
The latter leads to a limiting observable count rate for 
a given R.A. and decl. Using Xspec, we can compute 
the Galactic absorption-corrected flux limit versus sur- 
vey area for the RASS-AGN based on count rates, Nn 
values, and L = 2.5. We then compute the comoving 
volume (zi ow = 0.16) available to each object (V a ) for 
being included in the sample following lAvni fc Bahcalll 
(1980), usi ng the Galactic absorp tion-corrected object 
fluxes from lAnderson et al.l (|2007f ). From this we calcu- 
late the comoving number density as /jagn = J2i V^M; 
where i sums over each object. 

2.3. Defining a Common Survey Geometry of the 
RASS-AGN and LRG Samples 

The RASS/SDSS lAnderson etafl (|200l sample is 
based on the SDSS DR5, while the LRG sample is drawn 
from SDSS DR7. We make use of DR7 for the LRG 
sample as it contains the latest available and furthest 
advanced version of the SDSS products. Numerous cor- 
rection have been applied in comp arison to earlier data 
releases (see lAbazaiian et"aLll2009l ; e.g., updated photo- 
z, repeated observations for few regions with poor seeing 
in previous data releases, filling holes in DR6 region, cor- 
rection of instability in the spectroscopic flat- fields) . We 
then limit our LRG sa mple to the region covered by the 
lAnderson et al.l (|2007h AGN sample for the CCF calcu- 
lation. 

The SDSS survey geometry and completeness 
are expressed in t erms of spherical polygons 
([Hamilton fc Tegmarkl l2004f ). Publicly available ge- 
ometry and completeness files are not available for DR5, 
which would have the largest common survey area be- 
tween the LRG and the RASS-AGN samples. Therefore, 
we use the latest version available prior to DR5: the 
DR4+ geometry file, 5 , whic h is a subset dPark et al 
2001 of the SDSS DR5 (jAdelman-McCarthv et al 
2007ft and covers 5540 deg 2 (DR5: 5740 deg 2 ) 



The final, fiber collision-corrected LRG sample used 
here is based on DR7 but reconfigured to include only 
the DR4+ survey area that has DR7 completeness ra- 
tios of /compi > 0.8. The corresponding area of this 
sa mple is 5468 deg 2 . T his reconfiguration of the to- 
tal lAndersoneOLrilQQS) sample from DR5 to DR4+ 
eliminates 287 broad emission line objects, leaving 5937 
AGNs. Applying the redshift range selection of the LRG 
sample results in the number of objects given in Table [1] 

5 http: //sdss .physics .nyu. edu/lss/dr4plus 



for each AGN sample. Figure [2] shows the sky coverage 
of our final RASS-AGN sample and the LRG sample re- 
configured to a common DR4+ geometry which is used 
for the calculation of the CCF. 

3. MEASURING THE CROSS-CORRELATION FUNCTION 

A commonly used technique for measuring the spatial 
clustering of a class of ob ject s is the two-point correlation 
function £(r) ([Peebles! 1 1980lh which measures the excess 
probability dP above a Poisson distribution of finding 
an object in a volume element dV at a distance r from 
another randomly chosen object: 



dP = n[l + £(r)]dV, 



(2) 



where n is the mean number density of objects. The ACF 
measures the excess probability of finding two objects 
from the same sample in a given volume element, while 
the CCF measures the excess probability finding an ob- 
ject from one sample at a distance r from another object 
drawn from a different sample. The two-point correlation 
function, £(r), is equal to for randomly distributed ob- 
jects, and £(r) > if objects are more strongly clustered 
than a randomly distributed sample. 

In practice, the correlation function is obtained by 
counting pairs of objects with a given separation and 
comparing to the number of pairs in a random sample 
for the same separation. Differe nt correlation estimators 
are described in the literature. iDavis fc Peebles! (|1983l ) 
give a simple estimator with the form 



DD(r) 
DR{r) 



- 1, 



(3) 



where DD(r) is the sum of the data-data pairs at the sep- 
aration r and DR(r) is the number dat a-random pairs; 
both p air counts have been normalized. lLandv fc Szalayl 
(|1993f) suggest a more advanced estimator: 



C= — [DD-2DR + RR], 



(4) 



where RR is the normalized number of random-random 
pairs; this estimator yields errors similar to what is ex- 
pected for Poisson errors only. 

Because we measure line-of-sight distances from red- 
shifts, the measurement of £ is subject to the redshift- 
space distortions due to peculiar velocities. To separate 
the effects of redshift distortions, the spatial correlation 
function is measured as a function of two components 
of the separation vector between two objects, i.e., one 
perpendicular to (r p ) and the other along (n) the line of 
sight. Therefore, £(rp, it) is extracted by counting pairs 
on a 2D grid of separations r p and n. The real-space 
correlation function £(r) can be recovered by integrat- 
ing along the 7r directi on and computing the projected 
correlation function by IDavis fc Peebles! (|1983[ ) 



w p (r p ) = 2 / d7r£(r p ,7r) 



= 2 



rf^[(^ + y 2 ) (1/2) 

\drZ(r) (r 2 -rlY 1/2 



(5) 



which is independent of redshift-space distortions. The 
variable y represents the real-space separation along the 
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Fig. 2. — Spatial distribution of the total RASS-AGN sample (black squares) and SDSS LRG sample used here (dots, our LRG sample) 
restricted to the SDSS DR4+ geometry. The total covered area is 5468 deg 2 . 



Table 1. Properties of the LRG and RASS-AGN Samples. 



Sample 
Name 


2-range 


M|= - 3 Range 
(mag) 


Number 


< "LRG > 

(h 3 Mpc" 3 ) 


(*) 


< M 9 Z=0 - 3 > 
(mag) 


LRG sample 


0.16 < z < 0.36 


-23.2 < M|= ' 3 < -21.2 


45899 


9.6 x 10~ 5 


0.28 


-21.71 




2-range 


r 0.1-2.4keV 
X 

Range (erg s ) 


Number 


< "AGN > 

(h 3 Mpc" 3 ) 


<*> 


, r 0.1-2.4kcV . 
< ^ x > 
(erg s 1 ) 


Total RASS-AGN sample 
Low L x RASS-AGN sample 
High L x RASS-AGN sample 


0.16 < z < 0.36 
0.16 < 2 < 0.36 
0.16 < 2 < 0.36 


L x < 1-95 x 10 44 
L x > 1.95 x 10 44 


1552 
990 
562 


6.0 x 10- a 
5.8 x 10~ 5 
1.2 x 10~ 6 


0.25 
0.24 
0.28 


1.49 X 10 44 
8.81 x 10 43 
3.78 x 10 44 



line of sight. For a power law correlation function 

£(r)= (^) \ (6) 

ro and 7 are readily extracted from the projected corre- 
lation function using the analytical solution 

,. (r) _ r (roY r(l/2)r(( 7 -l)/2) 

wA r p ) r p y r J Th/2) , (7) 

where T(x) is the Gamma function. 

The aim of this paper is to more accurately mea- 
sure the clustering properties of low- z AGNs than has 



been measured previously using ACFs fe.g.. lMullis et al.l 
120041 : iGrazian et al.ll2004T ): we accomplish this by mea- 
suring the CCF of the AGNs with higher- density LRGs 
in the same volu me. Assuming a linear bias, we follow 
iCoil et al.l (|2009l) and infer the ACF of the AGN sample 
using 



w p (AGN|AGN) 



w p (AGN|LRG)] s 
w p (LRG|LRG) 



(8) 



where w p (AGN|AGN) and w p (LRG|LRG) are the ACFs 
of the RASS-AGNs and the LRGs, respectively, and 
w p (AGN I LRG) is the CCF of the RASS-AGNs with 
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the LRGs. The LR G ACF is studied extensively in 
IZehavi et al.l (|2005a| L where the estimator in Equa- 
tion Q is used; the results are given here in Table [21 
The RASS-AGN-LRG CCF is computed here using the 
estimator given in Equation ([3]) 



-Dagn ^lrg 



AGN-LRG 



- 1. 



(9) 



-Dagn -Rlrg 

We use this estimator as it requires a random catalog 
for only the LRG sample (-Rlrg), which is homoge- 
nous, volume-limited, and has a well-understood selec- 
tion function. The estimator given in Equation (|4]) would 
require a random catalog for the RASS-AGN sample, 
which would be subject to possible systematic biases 
due to difficulties in accurately modeling the position- 
dependent sensitivity limit. Especially, the changing 
Galactic absorption over the sky causes variations in the 
flux limit, which would require spectrum-dependent cor- 
rections. 

3.1. Construction of the Random LRG Sample 

The generation of random samples is crucial for a 
proper measurements of the correlation function. The 
objective is to construct a sample of randomly dis- 
tributed sources that have the same observational survey 
biases as the real sample. Use of the estimator given in 
Equation (j3|) requires that the construction of a random 
sample for only the LRG population is needed. 

The LRG sample used here has been corrected for fiber 
collisions (Section I2.1.2j) : therefore, we do not have to 
consider this bias in the construction of the random LRG 
sample. The SDSS survey geometry and completeness 
ratio for a given field are given by the SDSS geometry 
files. For a set of random R.A. and decl. values, they al- 
low us to determine if an object is covered by the SDSS 
DR4+ geometry and the spectroscopic completeness ra- 
tio at that location. Only objects covered in DR4+ with 
a DR7 spectroscopic completeness ratio / CO mpi > 0.8 are 
accepted for the random sample required here. Addition- 
ally, the DR7 spectroscopic completeness is used as the 
probability that an object is kept for the random sam- 
ple. If the completeness ratio is 0.9, the object has a 
90% chance of being included in the final random LRG 
sample. This procedure takes into account the fact that 
survey regions with a high spectroscopic completeness 
ratio have, on average, a higher object density than less 
complete areas. 

The corresponding redshift for a random object is as- 
signed based on the smoothed redshift distribution of the 
LRG data sample. Th e smoothing has been ma de by ap- 
plying a least-squares (jSavitzkv fc Golayill964l ) low-pass 
filter to the observed LRG redshift distribution. We com- 
pare the redshift distribution of the LRG data sample, 
its smoothed profile, and the redshift distribution of the 
random LRG sample in Figure [H The shape of the LRG 
redshift distribution is caused by the superposition of two 
selection criteria for the LRGs. At low redshifts, most 
of the LRGs are selected in the SDSS main galaxy sam- 
ple. This flux-limited selection reaches the maximum at 
z ~ 0.22. The 'cut F LRG selection provides objects 
already from the lowest redshifts with a fast increasing 
number of objects for higher redshifts. This selection 
reaches its flux limit at z ~ 0.35. The 'cut IF selection 
plays an important role only at z > 0.42. The selection 




0.25 
redshift 

Fig. 3. — Redshift histogram of the LRG data sample (solid black 
lines), the smoothed profile (red), and the random LRG sample 
(blue-dashed, renormalized to the total number of LRGs). 



dependence in diff erent redshifts is illust rated in detail 
in Figure 12-14 of lEisenstein et al.l (|200lD . 

The random catalog contains 100 times as many ob- 
jects as the LRG sample. This value is chosen to have 
an adequate number of pairs in the -Dagn -Rlrg sample 
at the smallest scales measured here. 

3.2. Errors and Covariance Matrices 

The calculation of realistic error bars on measurements 
of the correlation function has been a subject of debate 
since the earlie st measuremen t s. Di fferent methods are 
summarized in iNorberg et a l. (2009). Adjacent bins in 
w p (r p ) are correlated, as are their errors. The construc- 
tion of a covariance matrix Mij , which reflects the degree 
to which bin i is correlated with bin j, is needed to obtain 
meaningful power law fits to w p {r p ). 

We estimate the statistical errors of our correlation 
measurements using the jackknife method. We divide the 
survey area into A^t = 100 sections, each of which is ^55 
deg 2 . We calculate w p (r p ) Nt times, where each jack- 
knife sample excludes one section. These At jackknifc- 
resampled RASS-AGN ACFs are used to derive the co- 
variance matrix Mj» by 



Mi. 



N T -1 



E 

Lfc=l 



w k (r p ,i)- < w{r Pti ) > 



x\Wk(r P ,j)- < w(r P j) > 



(10) 



where Wk{r p ^) and uik{r p j) are from the fc-th jack- 
knife samples of the RASS-AGN ACF and < w(r P}i ) >, 
< w(r P: j) > are the averages over all of the jackknife 
samples. The la error of each bin is the square root of 
the diagonal component of this matrix (cr^ = \J Ma ) . To 
calculate the covariance matrix of the RASS-AGN ACF, 
which is determined using Equation ([8]) , we compute the 
RASS-AGN ACF for each of the At jackknife samples 
from the corresponding RASS-AGN-LRG CCF and LRG 
ACF of each jackknife sample. 

3.3. The RASS-AGN Auto-correlation Function 
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Table 2. Results of Power Law Fits to the LRG ACFs and 
AGN-LRG CCFs 



Sample 


r (h,- 1 Mpc) 


7 


Our LRG sample 
Zehavi subsamplel 


9.80±0.20 


1 qfi+002 
1.94±0.02 


Total RASS-AGN sample 
Low L x RASS-AGN sample 
High L x RASS-AGN sample 


6 93+^'' 

6 12+ - 50 
°' lz -0.53 

7 74 +0.40 
' ■ ' 4 -0.43 


4-8O-0.04 
1 94+0.10 
1,3 -0.08 
i OO+0-08 
l- yz -0.08 



Note. - Values of r$ and 7 are obtained from a 
power law fit to w p (r p ) over the range r p =0.3-40 h 
Mpc for all samples using the full error coy ariance matrix. 
For the LRG ACFs, IZehavi et ail (|2005al) and we used a 
7r max = 80 h- 1 Mpc, while for all CCFs 7r max = 40 h' 1 
Mpc was applied. 



We compute the CCF between RASS-AGNs and LRGs 
for the total sample, the low Lx sample, and the high Lx 
sample (see Table 2), as well as the ACF of the LRGs. 
We measure rp in a range of 0.3-40 h^ 1 Mpc in 11 bins 
in a logarithmic scale, while tt is computed in steps of 
5 h- 1 Mpc in a range of tt = - 200 h- 1 Mpc. The 
resulting £(rp, tt) are shown in Figured] for r p = — 30 
h^ 1 Mpc and tt = - 80 h^ 1 Mpc. Note the flattened 
contour at tt ~ 40 hr 1 Mpc in the LRG ACF. This is the 
first direct observation of the c oherent infal l for LRGs as 
expected by the Kaiser effect (|Kaiserlll987j ). 

Although Equation ([5]) requires an integration over tt 
to infinity, in practice an upper bound of integration 
(TTmax) is used: 

w p (r p ) = 2 / dTT^(r p ,Tr) . (11) 
Jo 

The value of 7r max has to be large enough to include most 
correlated pairs and give a stable solution, but not be so 
large as to unnecessarily increase the noise in the mea- 
surement. To determine the appropriated 7r max values 
for our correlation functions, we determined the correla- 
tion length ro for a set of 7r max values by fitting w p {r p ) 
with a fixed 7 = 1.9 over a r p range of 0.3-40 hr 1 Mpc. 

Figure [5] shows that the LRG ACF saturates at 7r max = 
80 hr 1 Mpc. The changes in the correlation lengths 
above this value are wel l within the uncertainties. There- 
fore, as in lZehavi et al.l (|2005a[) . we use an upper bound 
of the LRG ACF integration of 7r max = 80 hr 1 Mpc. Ta- 
ble [5] shows the values of ro and 7 f or a power law fit in a 
range of 0.3-40 h' 1 Mpc (as used in lZehavi et al.ll2005aD . 
Both results well agree within their uncertainties. 

The CCFs between the different RASS-AGN samples 
and the LRGs saturate at 7r max = 40 hr 1 Mpc (Fig- 
ure [5]). At higher values of 7r max , the signal-to-noise ratio 
degrades and no significant change in ro occurs. There- 
fore, we use 7r max = 40 h^ 1 Mpc as an upper bound of 
integration for all AGN-LRG CCFs. The difference in 
the measured LRG ACF using 7r max = 40 hr 1 Mpc and 
TTmax = 80 h~ 1 Mpc is only 3%. We expect that the 
growth of the CCF between 7r max = 40 hr 1 Mpc and 
TTmax = 80 h^ 1 Mpc is about the same order. Since this 
is much smaller than the errors in the CCF, a use of 
TTmax = 40 hr 1 Mpc for the CCF is reasonable. The cor- 
relation length of the CCFs of the different RASS-AGN 



samples and the LRGs is given in Table [2j 

As the area covered by SDSS DR4+ is not contigu- 
ous (see Figure [5]) , we computed the CCF for different 
subsamples of the SDSS DR4+, to check that there were 
no biases introduced by using non-contiguous regions of 
the sky. Excluding survey areas with R.A. < 70 ° and 
R.A. > 300°, the isolated area around R.A. ~ 260° and 
deck ~ 60°, the somewhat patchy areas at 190 ° < 
R.A. < 250° and 0° < deck < 35°, and combinations 
thereof, results in measurements of the CCFs that all 
agree well within their uncertainties. Changing the step 
size of tt from 5 to 2.5 h^ 1 Mpc alters w p (r p ) by a neg- 
ligible amount. 

Instead of using the derived values of ro from the power 
law fits of the LRG ACF and AGN-LRG CCFs to com- 
pute the RASS-AGN ACF, we use Equation ([8} and use 
the full w p (r p ) functions. Figure |5] shows w p (r p ) for the 
RASS-AGN ACF, the LRG ACF, and the RASS- AGN- 
LRG CCF. FigureEshows the CCFs for the low and high 
Lx RASS-AGN samples. 

We fit power laws to the ACFs of the different RASS- 
AGN samples. The fit uses the covariance matrix and 
minimizes the correlated \ 2 values according to 

JVbins N bina , v 

x 2 =EE K(^,o-< ode W)) 

i=l j=l ^ ' 

XM^L^)-^'^)) (12) 

We only fit the data in a range r p =0.3-15 h~ x Mpc, as 
the clustering signal above r p =15 Mpc is not well- 
constrained for the low Lx RASS-AGN sample. The up- 
per end of r p has also been chosen because we will later 
convert the fit results into cts.agn which involves only the 
pairs within 16 hr 1 Mpc. Contour plots of the resulting 
values of ro and 7 are shown in Figure [8] for the different 
RASS-AGN samples. The derived best-fit values, as well 
as the best-fit ro values with a fixed power law slope of 
7 = 1.9, are given in Table [3] Based on the error on 
ro for a fixed 7 = 1.9, we estimate the clustering signal 
to be detected at a ~11ct, ~5<t, and ~8er level for the 
total, the low Lx, and the high Lx RASS-AGN sample, 
respectively. The difference in the clustering signal be- 
tween the low Lx, and the high Lx RASS-AGN sample 
is detected at the ~2.5cr level. 

The clustering strength is commonly expressed in 
terms of the rms fluctuation of the density distribution 
over the sphere with a comoving radius of 8 hr 1 Mpc. 
For a power l aw correlation func tion, this value can be 
calcu lated by (jMivaji et alJl2007t Section 59 of iPeeblesI 
[19801) 

Ka-) 2 = L 2 (7)( ¥ ^) 7 , (13) 

where J 2 {l) = 72/ [(3 - 7) (4 - 7) (6 - 7)2^]. The un- 
certainty of o"8,agn is computed from the ro versus 7 
confidence contour of the one-parameter fit based on a 
correlated % 2 = x min + 1.0 (Figured]). 

Based on erg, agn, we further calculate the RASS-AGN 
bias parameter b — (Js,AGN(z)/ag,(z). This quantity al- 
lows us to compare the observed AGN clustering to the 
underlying mass distribution from linear growth theory 
dHamTlton 200l]). We use a normalization to a value of 
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Fig. 4. — Contour plot of the CCF between RASS-AGNs and LRGs for the total sample, the low L x AGN sample, the high L x AGN 
sample, as well as the LRG ACF (left to right). Contour lines show constant correlation strength for the 2D correlation function t;(rp,ir). 
The data are not smoothed. The contour levels are 0.0 (dotted line), 0.1, 0.2, 0.5, 1.0 (thick solid line), 2.0, and 5.0. 



Table 3. Power Law Fits to the RASS-AGN ACF. 



Sample 
Name 


ro 

(h,- 1 Mpc) 


7 


^c,7=1.9 

(h- 1 Mpc) 


°"8,AGn(^) 


b(z) 

Kagn(*)M(«)) 


log M DM H 

(ft" 1 M Q ) 


Total RASS-AGN sample 
Low L x RASS-AGN sample 
High L x RASS-AGN sample 


, 9 o+0.44 
4.2S_ Q 54 

o o 9 +0.64 
°- oz -0.83 
5 44+0-71 


1 67+ - 13 
1 - ot -0.12 

l 7 O+0.40 

1 - ' J -0.27 
, oc+0.20 
l.Sb_ 21 


4 " 3Z -0.41 
q Qfj+0.58 
°- zo -0.69 
c .0+0.64 
0.71 


q 77 +0.07 
u - ' '-0.08 

O.K.™ 
0.98±°;1| 


i-u±S:i§ 
1-4418:5? 


19 co+0.20 
lz.jo_Q.33 

11.83+° X) 55 

13 10+ 024 
1J.1U_ 43 



Note. — Values are obtained from a fit to w p (r p ) in a range of 0.3-15 h 1 Mpc for all samples using the full error 
covariance matrix. 



f7g(z = 0) = 0.8 which is consistent with the Wilkin- 
son Microwave Anisotropy Probe (WMAP) DR5 for a 
ACDM+SZ+LENS model 6 . The errors on b are derived 
from the standard deviation of cts.agn- 

Because we measure the CCF to infer the ACF, the 
resulting effective redshift distribution for the cluster- 
ing signal is determined by both the redshift distribu- 
tion of the LRG sample and the RASS-AGN sample: 
Nccf(z) = N WG (z) * Arass-agnW- The median red- 
shift of N CC f(z) is zccf = 0.27, 0.24, 0.31 for the to- 
tal, the low Lx, and the high Lx RASS-AGN samples, 
respectively. The difference compared to the mean red- 
shift of the different RASS-AGN samples (see Table P 
is at most Az = 0.03. Our measurements of cts.agn and 
b{z) (using zccf) for all RASS-AGN samples are listed 
in Table _0 



// RASS-AGNs are hosted by typical L* galaxies, the 
fraction of L* galaxies hosting RASS-AGNs can be cal- 
culated using the observed number density of bot h host 
galaxies and RASS-AGNs. iBlanton et ait (|2003bf ) mea- 
sure the number density of L* galaxies at z = 0.1 to 
be <p*_ haad = 0.0149 ± 0.0004 h 3 Mpc" 3 . Using the ob- 
served number density in the total RASS-AGN sample 
(Table [J), we find that only ~0.4% of L* galaxies could 
harbor RASS-AGNs. 

3.4. Host Dark Matter Halo Mass 

Our clustering results can be used to estimate the "typ- 
ical" dark matter halo mass Mdmh hosting our different 
AGN samples. The dark matter halo mass is reflected 
by the bias parameter 6, which reflects the clustering 



amplitude relative to the u nderlying dark mat ter distri- 

3 |http : //lambda . gsf c ■ nasa ■ gov/product /niap/dr3/params/lcdm_s2_lfel^tiiaHp5Uaffil g Equation 8 of lSheth et alj ()200lD . we Com- 
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Fig. 5. — Correlation length ro vs. 7r max using Equation (0 
for the LRG ACF (triangles), the high L x RASS-AGN-LRG CCF 
(diamonds), the total RASS-AGN-LRG CCF (boxes; for illustra- 
tion purposes, this has been shifted +2.0 in 7r max direction), the 
low L x RASS-AGN-LRG CCF (crosses; shifted by +4.0 in vr max 
direction). 
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Fig. 8. — Probability contours for the power law normalization 
ro and slope 7 of the total RASS-AGN sample (black, central con- 
tour), low Lx RASS-AGN sample (gray, left contour), and high 
Lx RASS-AGN sample (gray, right contour). The solid contours 
represent the lcr confidence intervals (68.3%) for a one-parameter 
fit based on a correlated \ 2 = Xmin + 1-0 (Equation Q12I ). while 
the dotted lines illustrate the corresponding intervals for a two- 
parameter fit (lcr, correlated \ 2 = Xmin + 2.3). 
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Fig. 6.— Projected LRG ACF (boxes), RASS-AGN-LRG CCF 
(triangles), and the inferred RASS-AGN ACF (diamonds). For 
all data points, we show the lcr uncertainties. The obtained best 
power law fit of the RASS-AGN ACF, using the covariance matrix 
and fitting over r p = 0.3 — 15 h~ 1 Mpc, is shown as a dashed line. 
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Fig. 7. — CCFs of the low (stars) and high (diamonds) L x RASS- 
AGN samples are shown with their corresponding uncertainties. 



pute the expected large-scale Eulerian bias factor for dif- 
ferent dark matter halo masses. The required ratio of the 
critical overdensity to the rms fluctuation on a given size 
and m ass is calculated by v — 5 CT /a{M,z) (jSheth et alJ 
|2001| ). Assuming S CT ~ 1.69, we comput e o~(M, z) using 
Equat ions (A8), (A9), and (A10) from Ivan den Bosch 

20021) ; this approach was also used in iHickox et al. 

2009). The typical dark matter halo mass for the total 
RASS-AGN and the high L x RASS-AGN sample is found 
to be log (Mdmh/^Mq)) = 12.5812:12, 13.10lg;J|, 
respectively. For the low Lx RASS-AGN sample (log 
(Mdmh/^ 1 M )) = 11.83+ ' 55 ) we are not able to con- 
strain a lower limit on the dark matter halo mass. This 
is because our lower limit corresp onds to b = 0.71 a nd 
the minimum b value derived from lSheth et al.l (|2001[ ) at 
z = 0.24 is 6 mi „ = 0.72 (at M DM h ~ 10 9 3 h' 1 M & ). 

4. DISCUSSION 

4.1. Comparison to Other X-ray Clustering 
Measurements 

Our findings can be directly compared to previous at- 
tempts to measure the spatial ACF (using spectroscopic 
rcdshifts) of low-;? RASS-b a sed A GNs at similar X-ray lu- 
minosities. iGrazian et al.l (|2004f ) found a redshift-space 



clustering length of sq = 8.64 



+2.00 h -i 

-2.08 



h Mpc in the 



ASIAGO-ESO/RASS QSO survey when fitting their cor- 
relation function with 7 fixed at 1.56. They did not pub- 
lish results for fits of ro and 7. A second low- 2: RASS- 
AGN clustering stud y with very s imila r AGN proper- 
ties was conducted bv lMullis et alJ (|2004j ) in the ROSAT 
North Ecliptic Pole survey. Their best-fit values are 
r n = 7.5+H h~ l Mpc and 7 = 1.85+H2. Our inferred 



4.281°^ 



r a = L0ii; 2 ti " Mpc and 7 — jl. ou _ 80 
ACF of the total RASS-AGN sample yield r 

h ^ Mpc and 7 = 1.67lg Both studies used Pois- 
son errors which underestimate systematic effects. Us- 
ing the largest sample of X-ray selected AGNs ever ap- 
plied to clustering measurements, and by measuring the 
cross-correlation with more numerous galaxies, we reach 
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significantly lower uncertainties compared to the previ- 
ous studies and detect a clustering signal at the ~llcr 
level. We also find a m uch smaller c orrela tion length. 
iGrazian et ail (|2004 and lMullis et all (|2004D derive their 
results from fitting the ACF at much larger separations 
(s, rp values) than we do. However, considering their 
larger uncertainties, their correlation lengths differ from 
ours by only \a-1a . The low number of X-ray selected 
AGNs used in both studies did not allow them to split 
their samples in bins of either luminosity or redshift. 

At higher redshift, data from more modern X-ray tele- 
scopes such as XMM-Newton and Chandra have been 
used to measure the clustering of X-ray selected AGNs 
in various deep survey fields (Table 2]). These surveys 
are significantly more sensitive than R OS A T surveys but 
have sampled smaller comoving volumes (due to smaller 
sky area coverage) and may therefore be affected by cos- 
mic variance. The median X-ray luminosity of these sur- 
veys is roughly 1 mag fainter than the X-ray luminosi- 
ties detected in ROSAT-ba,sed catalogs. Furthermore, 
XMM-Newton and Chandra are sensitive at both soft 
energy bands (as ROSAT: 0.1-2.4 keV) and hard energy 
bands (2-10 keV). This influences the sample selection; 
ROSAT samples are dominated by X-ray unabsorbed 
type I AGNs, while XMM-Newton and Chandra samples 
consist of a nearly equal mix of type I AGNs and ab- 
sorbed type II AGNs. A direct comparison of ROSAT 
and XMM-Newton/ Chandra AGN clustering measure- 
ments is therefore challenging. 

The inferred typical host dark matter mass for RASS- 
AGNs at z — 0.27 is consistent with that found at higher 
redshift in CDFN and AEGIS, though other X-ray AGN 
clustering studies find significant larger values at these 
redshifts (i.e., CDFS, CLASXS, and COSMOS; see Ta- 
ble H for details) . 

One possible explanation of the differences in the 
clustering signals between the low and high-z results 
could be the presence of a large fraction of type II 
AGN s in the XMM-Newton a nd Chandra samples. How- 
ever, iGilli et all ([2005, 2009) measured the correlation 
strengths separately for soft and hard AGN samples, as 
well as for broad emission line versus narrow emission 
line AGN samples, and did not find a significant differ- 
ence. However, larger samples with smaller u ncertainties 
are needed to definitively test this question. IGilli et all 
(|2005l I2009D explain the large values of r Q in the CDFS 
and COSMOS by the cosmic variance, i.e., very promi- 
nent redshift spikes at z ~ 0.7 and z ~ 0.36, respectively. 
Removing these overdensities in the total X-ray sample 
results in correlation lengths of r = 3.8l 2 ; 7 , 6.1±0.8 
h^ 1 Mpc, respectively, similar to the values found here 
at lower redshift. This result underscores the need to 
include cosmic variance errors when using samples that 
do not probe a large volume. Among measurements that 
correctly account for cosmic variance, cosmic evolution 
and differences in X-ray luminosity may be responsible 
for the different clustering results in the literature. We 
will discuss this point in detail in Section [4.31 

In Table [H we compare the main properties of the 
different clustering measurements of AGNs selected at 
different wavelengths (e.g., X-ray and optical) and star- 
forming and quiescent (i.e., blue and red) galaxies at 
similar redshifts. Table|4]aims to compare different stud- 
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Fig. 9. — Bias parameter &agn = &8,AGn(z) /[&8D(z)] as a func- 
tion of redshift for various X-ray AGN selected samples (see the 
text for details). Gray triangles represent luminosity-dependent 
subsamplcs of the main AGN sample used here. The gray filled el- 
lipses indicate op t ical AGN clusterin g measurements (left t o right: 
ICoil et ail j2007h : IRoss et al.1 p009h : IPorciani et alj pOOl O. The 
dotted lines show the e xpected b(z) of ty pical dark matter halo 
masses Mdmh based on Shcth ct al. (2001). The masses are given 
in log Mqmh m units of h~ 1 Mq. 



ies as uniformly as possible to detect general trends in 
the clustering of AGNs to z ~ 1. Therefore only stud- 
ies that are based on the real-space correlation function 
£(r) are included here. Many results, especially for op- 
tically selected AGNs, measure only the redshift-space 
correlation function, £(s), and attempt to model redshift- 
space distortions to infer £(r). However, £(s) is not well- 
approximated by a power law and the derived results 
depend therefore on the fitted s range. Additionally, sys- 
tematic uncertainties in the modeling make direct com- 
parisons to £(r) measurements extremely difficult. To fa- 
cilitate comparison between samples, we recalculate the 
bias, b{z), for each study given in Table |4] in a consistent 
manner. Values of o^agn/gal are com P u ted using con- 
fidence contours in r -j space (where x 2 = Xbestfit + ■*-) 
published in the literature. If no confidence contours are 
given, we use the best- fit tq, 7 values and propagate their 
error to derive the la error on o^agn/gal- ^ e excm de 
papers that report only b values but do not publish ro and 
7 values. Where possible, we report the median redshift 
of objects in each study that contribute to the clustering 
signal. We use z c s whenever available, otherwise we use 
the median redshift (z), or the mean redshift ((2)), in 
that order, to compute erg(z) using our normalization of 
a$(z = 0) = 0.8 throughout. 

Figure [5] shows the derived bias factors b of various X- 
ray selected AGN samples as a function of redshift. The 
results shown in this figure suggest that high redshift, low 
X-ray luminous AGNs tend to reside in higher mass dark 
matter halos (Mdmh ~ 10 13 h^ 1 Mq) than low redshift, 
high luminous AGNs. Broad-line X-ray and optical AGN 
measurements seem to show a lower bias values than the 
narrow-line and/or fainter X-ray AGN measurements. 

4.2. Comparison to AGN Clustering Measurements at 
Different Wavelengths 



12 



Krumpe et al. 



Table 4. Comparison of Published Real-space Clustering Measure- 
ments of X-ray Selected AGNs and Optically Selected AGNs, as well 
as Red and Blue Galaxies. 



Fitted Area Object 



L,M 



Band 



Bias 



Sample 


(h^ 1 Mpc) 


7 


Range 


(deg 2 ) 


Number 


(erg/sj,(mag) 


(keV), filter 


z 


b 


Ref. 


RASS 

Low-L RASS 
High-L RASS 


i.28+°f 4 

o o 9 +0.64 
— U.oo 

— u.ya 


, 7 o+0.40 

27 

i.86±g-;5° 


0.3-15 
0.3-15 
0.3-15 


5468 
5468 
5468 


1552 
990 
562 


1.4 x 10 44 
9.8 x 10 43 
3.4 x 10 44 


0.1-2.4 
0.1-2.4 
0.1-2.4 


0.27 
0.24 
0.31 


1.11+° 10 
1.44+S : f 7 


This work 
This work 
This work 








X-ray Selected AGN Clustering Measurements 








NEP 


7 5+ 2 ' 7 
1 ->J_4.2 


1 otr + 1.90 

±.Od_Q g Q 


5-60 p 


80.7 


219 


9.2 x 10 43 


0.5-2.0 


0.22 B 


t-»3_ 61 


Mu04 


CDFN 


5.5 ±0.6 


1.50 ±0.12 


0.2-10 p 


0.13 


160 


1 x 10 43 


0.5-10 


0.96 


-, o 7 +0.14 
l-8'_0 16 


Gi05 


CDFS 


10.3 ±1.7 


1.33 ± 0.14 


0.2-10 p 


0.1 


97 


1.6 x 10 43 


0.5-10 


0.84 


2 64+ ' 29 
z ' D4 -0.30 


Gi05 


CLASXS 




2.1 ± 0.5 


1-30 


0.4 


233 


9 x 10 43 


2-8 


1.2 


o tro + 2.49 
- °-1.38 


Ya06 


CDFN 


5.8+|-° 

— 1.5 


-, Qo+0.12 
— 0. 14 


0.2-15 


0.13 


252 


3 x 10 42 


2-8 


0.8 


1.77+°'?° 

—0.15 


Ya06 


AEGIS 


5.95 ± 0.90 


1.66 ± 0.22 


0.1-8 


0.63 


113 


6.3 x 10 42 


2-10 


0.90 


1.97 + °' 2 « 

—0.25 


Co09 


COSMOS 


—0.48 


1 »H+ 006 
—0.07 


0.3-40 p 


1.96 


538 


6.3 x 10 43 


0.5-10 


0.98 


3.08+°'H 

— 0. 14 


Gi09 








Optically Selected AGN Clustering Measurements 








2QZ 


4.8i; : g 


, co+0.19 
l. jo_q 21 


0.8-20 


~700 


13989 


-23.82 




1.47 s 


o n7+0.35 
Z - | J'_0.44 


Po04 


2QZ-subsl 


5.4i?-» 


r> n r> + 0.36 

z - uz _0 33 


2-20 


~700 


4928 


-23.13 


bj 


1.06 B 


2 14+°- 71 
—0 55 


Po04 


2QZ-subs2 


4.3±|;« 


1 4Q+032 


2-20 


~700 


4737 


-23.84 


bj 


1.51 B 


1 qi+0.73 

l.»0_Q gQ 


Po04 


2QZ-subs3 


7.6+_H 


-, 7q +0.25 
1 - ,a -0.29 


2-20 


~700 


4324 


-24.30 


b] 


1.89 B 


o 71 +0.97 
°' 1 1 — 1.03 


Po04 


DEEP2 


3.1 ±0.6 


1.8(fixed) 


0.1-10 


3 


52 


-23.0 


B 


0.99 


, 1s +0.20 
t-tO-O 20 


Co07 


SDSS 


°- 40 -0.45 


x - JU -0.03 


1-130 


~4000 


30239 


~-25.8 


i 


(1.27) 


2 26+°' 14 
z - zo -0.18 


Ro09 








Galaxy Clustering Measurements 






1 41 +0.08 
1 -^ ll -0.08 




red-2dF 


6.10 ±0.34 


1.95 ± 0.03 


0.2-20 


~700 


36318 


1.26 


L* 


0.11 


Ma03 


blue-2dF 


3.67 ±0.30 


1.60 ±0.04 


0.2-20 


~700 


60473 


0.95 


L* 


0.11 


n ocr + 0.05 
U - 8t, -0.05 


Ma03 


red-SDSS 


5.67 ±0.37 


2.08 ±0.05 


0.1-10 


2497 


5804 


(-19.5) 


r 


(0.05) 


1.4x1°:°! 


Ze05b 


bluc-SDSS 


3.63 ±0.16 


1.69 ±0.04 


0.1-10 


2497 


8419 


(-19.5) 


r 


(0.05) 




Ze05b 


LRG-SDSS 


9.80 ±0.20 


1.94 ±0.02 


0.3-30 


~3800 


29298 


(-21.63) 


9z=0.3 


(0.28) 


2 5fi+°-° 6 


Ze05a 


red-AGES 


5.3 ±0.2 


2.1 ± 0.1 


0.3-10 


7 


3146 


-21.3 


r z =o.l 


0.41 


1 KO+0-06 
i - oy -0.06 


Hi09 


bluc-AGES 


3.8 ±0.2 


1.6 ±0.1 


0.3-10 


7 


3116 


-21.0 


r z =o.i 


0.38 


1 n7 +0.05 
J - u ' -0.05 


Hi09 
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Note. — Values for the fitted range (Column 4) are in units of h~ x Mpc. Values listed for L and M are the median 
values whenever these quantities were available. Otherwise, mean values are given and denoted by angled brackets 
(()). Superscript "P" on the fitted range indicates that only Poisson errors have been used for the fit, which causes 
a significant underestimation of the uncertainties. All other studies used either the bootstrap or jackknifc method to 
estimate their uncertainties. Superscript shows that the effectiv e redshift was g i ven in the stu dy, as oppose d to th e 
median redshift. Abbreviation: subs-subsample . Refer ences: G r04 Grazian et al. (2004); Mu04- jMullis ct al. ( 2004T) : 



GiOfilGilli et all (l2005h: Y aOB lYang et al.l (f200fih: Co09lCoil et all (f2009h: GiOQlGilli et al.l (l2009h: HiOQlHickox et al.l 
(120091) : Po04lPorciani et all (12004ft : Co074Coil et al.l (120071 ): Ro09 lRoss et alJ d2009f ): Ma03 lMadgwick et al l (I2003D: 



Ze05b 4zehavi et all (l2005b¥ Ze05a 4zehavi et alJ (|2005aE ~Co08- [Coil et all (|200cP leb6 4Meneux et al.l (|2006). 



The availability of large optical data sets such as SDSS, 
2dF QSO Redshift Survey (2QZ), and 2dF-SDSS LRG 
and QSO (2SLAQ) led to a series of large-scale clustering 
studies. However, many authors report only £(s) (e.g., 

ICroom et alll2005t Ida Angela et al.ll2008t ). which makes 
it difficult to compare their results to our £(r)-based re- 
sults. Other studies compute t he CCF between opticall y 
selected QSOs and LRGs fe.g.. lMountrichas et aT]|2009D . 
but do not infer the QSO ACF. Here we compare only 
to studies that measured £(r), listed in Table [4] (mid- 
dle section) . Most of the optically selected AGN samples 
contain samples that are orders of magnitu de larger than 
X-ray selected samples. ICoil et all (|2007l ) use the CCF 
between ~5000 DEEP2 galaxies and 52 QSOs to infer 
the QSO ACF. 



AGNs selected at optical wavelengths have relatively 
low clustering scale lengths, found over a wide range 
of redshifts. This translates into low bias parameters, 
not much larger than what we find for our low-z RASS- 
AGNs (Figure [10]). Radio- loud AGNs are found to 
be significa ntly more clustered than optically selected 
AGNs (e.g..lMagliocchetti et al.11200 1: Iffickox et al.ll2009l: 
iMandelbaum et al.l l2009f K Optically brig ht and radio- 
loud quasars are also luminous in X-rays (j Wilkes et al.1 
19941) and reside in very massive dark matter halos (e.g., 
Shen et al.ll2009h . 

Comparing the general trends seen in Figure [9] and 
[TOl X-ray and optically selected AGNs show similar clus- 
tering strength at z < 1. Optically selected AGNs at 
higher redshifts are found to have a lower clustering 
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Fig. 10. — Similar to Figure[9] now comparing our results to AGN 
samples selected at optical and radio wavelengths. Our results are 
shown at z = 0.3; gray triangles represent subsamples of the main 
AGN sam ple used here; gray b oxes show subsamples of the AGN 
sample by .Porciani et al. (2004). 
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Fig. 11. — Similar to Figure [9] now comparing our results to 
red and blue galaxy samples. For all galaxy results, the data point 
with lower b value corresponds to blue, star-forming galaxies while 
red galaxies have higher b values. 



than X -ray selected AGNs. The sample of lPorciani et al.l 
(|2004h has a higher mean lu minosity and highe r cluster- 
ing signal than the sample of iCoil et al.l (|2007|) . Within 
the errors, the clustering of optically selected AGNs is 
known to be in accordance with a redshift-independent 
dark matter halo mass o f M DM h ~ 10 12 — 10 13 h~ x 
Ma, (jPorciani et al.l [200l ICoil et al.l 120071 iRoss et al.l 
[200l . Low-luminosity AGNs (-23 < M 0Ar < -17) 
in the local universe have a typica l dark matter halo 
mass of M D mh ~8x 10 11 h^ 1 M Q (jMandelbaum et ail 
2009) , which is also c onsis tent with the M dmh given by 



Porciani et ail (|2004j ) and lCoil et aLll2007l for the high- 



redshift luminous AGNs considering the uncertainties in 
their clustering measurements. 

4.3. Luminosity Dependence of the Clustering Signal 

In the hierarchical model of structure formation, more 
massive galaxies should reside in more massive dark 
matter halos and therefore be more strongly clustered. 
More massive galaxies are expected to be more lumi- 



nous, which should lead to a luminosity dependence of 
the clustering signal. This result has been confirmed by 
clusterin g measurements of, e.g., SDSS galaxies at low 
redshifts l Zehavi et al.ll2005bl) and of, e.g., DEEP2 galax- 
ies at z = 1 (jCoil et al.ll2006l ). Whether this relation ap- 
plies also to the AGN luminosity is unclear. The AGN 
luminosity depends on the SMBH mass, mass accretion 
rate, and radiative accretion efficiency. In the simple 
case in which all SMBHs have the same Eddington ra- 
tio and the same dependence of radiative efficiency on 
the Eddington ratio, then higher-luminosity AGNs will 
have higher SMBH masses. // there is a correlation be- 
tween SMBH mass and dark matter halo mass, then the 
higher-luminosity AGNs will also be more strongly clus- 
tered. The observed correlation between the mass of the 
SMBH and t he stellar velocity dispersion in the bulge 
of the galaxy (jGebhardt et al.ll2000t iFerrarese fc Merrittl 
2000) suggests that a luminosity (X-ray and optical) clus- 
tering dependence could be feasible and physically moti- 
vated, since massive bulges preferentially reside in mas- 
sive DMHs. Alternatively, even among galaxies with the 
same SMBH mass, those that reside in a denser large- 
scale environment may have a higher chance of large mass 
accretion through galaxy mergers and interactions. How- 
ever, this trend has not been detected, possibly due in 
part to the relatively large uncertainties in the cluster- 
ing measurements of X-ray and optically selected AGN 
samples. 

Several studies have attempted to measu re the de- 
pende nce of clustering on X-ray luminosity. iGilli et all 
(2009) compute the clustering signal for a high and low 
X-ray luminosity AGN sample in the XMM-Newton- 
COSMOS field at a dividing line of j^-io faV = 1Q 44 
ergs - 1 and find no significant difference. lYang et"aL1 
(2006) studied an X-ray luminosity dependence in the 
Chandra se l ected AG N samples. As n oted by both 
lYang et al.l (|2006l ) and IGilli et"all (|2009ft . splitting the 
samples in different X-ray luminosities leads to the study 
of a high and low redshift sample of X-ray selected AGNs. 
These studies of the dependence of clustering on lumi- 
nosity are therefore hampered by the redshift-luminosity 
degeneracy. Taking into account the possib l e reds hift 
evolution of the clustering signal, lYang et aLI (|2006l) do 
not find a signi ficant dependence of clustering on the X- 
ray luminosity. ICoil et a l. (2009) measure the clustering 



of AGNs with 

-10kcV 



10 42 ergs 

< 10 42 ergs -1 at similar redshift (z = 0.8) 
and find no significant difference in their clustering prop- 
erties as well. 

Our sample here covers a limited redshift range, 0.16 < 
z < 0.36, such that we do not expect a significant contri- 
bution of the redshift evolution to the clustering signal 
within our sample. In addition, our use of the CCF with 
a large tracer set of LRGs yields an AGN ACF with high 
precision. Splitting the sample into a low and high X- 
ray luminosity sample at the commonly used AGN/QSO 
dividing line of i£ 5 ~ 10kcV = 10 44 ergs" 1 allows us to de- 
tect an X-ray luminosity dependence of the clustering at 
the ~2.5<7 level for the first time, in that X-ray luminous 
RASS-AGNs at low- z cluster more strongly than their 
low luminosity counterparts. 

Although, as mentioned above, the tracer LRG sample 
is approximately volume-limited, in order to clearly test 
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whether the difference in the clustering signal between 
the high and the low Lx RASS-AGN sample could be 
caused by a change of the LRG ACF with redshift (which 
could mimic a luminosity dependence of the inferred 
RASS-AGN ACF), the redshift de pendence of the LRG 
sample ACF must be analyzed. IZehavi et al.l ([2005a) 
study dependence of the LRG clustering amplitude with 
luminosity and redshift and split their 0.16 < z < 0.36 
LRG sample into 0.16 < z < 0.23 and 0.23 < z < 0.36 
subsamples; they find no redshift evolution in the LRG 
clustering signal. We repeat this test and split our LRG 
sample into two subsamples containing object in redshift 
bins of 0.16 < z < 0.26 and 0.26 < z < 0.36. We verify 
that the LRG ACF of both subsamples agree to much 
better than la and therefore exclude redshift evolution 
of the LRG clustering signal as a possible source of the 
detected luminosity dependence of the RASS-AGN clus- 
tering. 

Optical QSO/AGN samples based on SDSS, 2QZ, and 
2SLAQ samples yield comparable low uncertainties in the 
measured ACF due to the large sample sizes (up to a few 
10,000 objects), allowing measurements of the luminosity 
dependen ce of clustering in optically selected QSO/AGN 
samples. iMountrichas et ail (|2009T ) compute the CCF 
at z < 1 between QSOs and LRGs from 2SLAQ, 2QZ, 
and SDSS to break the redshift-luminosity degeneracy. 
They find little QSO-LRG cross-clustering dependence 
on QSO luminosity, implying dark matter halo masses 
of Mdmh ~ 10 13 h~ x Mq approximately independent of 
QSO luminosity. Da Angela et al. (2008) confirm that 
QSOs of different luminosities reside in dark matter halos 
of similar mass (Mdmh ~3x 10 12 hr 1 Mq) for z ~ 1.5. 

Using the optical luminosity in Table 2] and the a ox 
relation between the optical rest-frame flux density at 
2500 A and the X-ray luminosity flux density at 2 keV 
measured for 6224 AG Ns with broad emission lines from 
lAnderson etal] (j2007ft . we derive the median X-ray lu- 
minosities of the optical AGN samples. This allows us 
to compare these results with X-ray selected AGN clus- 
tering measurements. The optically selected AGN sam- 
ples have median X-ray luminosities of x^ 1_2 4kcV ~ 
3 — 9 x 10 44 ergs" 1 , which is very similar to the X-ray 
luminosity of our high Lx RASS-AGN sample. These 
samples are shown in Figure|9]as gray filled ellipses. Con- 
trary to our findings for the X-ray luminosity dependence 
at low-z, at higher redshift high luminosity optically se- 
lected AGNs appear to be less clustered than low lu- 
minosity X-ray selected AGNs. This implies that high 
luminosity AGNs at these redshifts lie in less massive 
dark matter halos than low luminosity AGNs. However, 
Chandra and XMM-Newton X-ray selected AGNs (which 
define the X-ray samples at higher redshift) contain a sig- 
nificant fraction of type II AGNs, while optically selected 
AGN samples contain mainly type I AGNs with broad 
optical emission lines. Therefore, X-ray and optically se- 
lected samples at high redshifts probe different kinds of 
AGNs. Optically selected, high luminosity AGN samples 
at high redshifts and high Lx RASS-AGN samples at low 
redshifts likely contain the same type of AGNs because of 
both ROSATs soft energy selection being biased toward 
X-ray unabsorbed type I AGNs and the similar X-ray lu- 
minosities measured for both samples. In terms of their 
clustering properties, optically selected AGN samples at 



higher redshift reside in host dark matter halo masses 
between M DM h ~ 10 12 h' 1 M Q and M DM h ~ 10 13 h' 1 
Mq, similar to the full range observed for the low red- 
shift RASS-AGN samples. Stated differently, the dark 
matter halo masses of both the low Lx and high Lx 
RASS-AGN samples are consistent with those of the op- 
tically selected, high luminosity AGNs at higher redshift 
(which have relatively large clustering uncertainties). 

To summarize the clustering trends seen, at low red- 
shift within the RASS-AGN broad-line sample, higher lu- 
minosity AGNs cluster more strongly than their lower lu- 
minosity counterparts, while at higher redshifts (z > 0.7) 
the higher luminosity optically selected AGNs are less 
clustered than the lower luminosity X-ray selected AGNs. 
However, the higher redshift X-ray AGN samples contain 
a large fraction of narrow-line AGNs, and so these trends 
are not necessarily comparing the same kinds of AGNs at 
all redshifts and luminosities. For example, the low Lx 
AGNs studied at high redshifts have lower luminosities 
than the low redshift, low L x RASS-AGNs (see Table gj. 
We propose two possible explanations for the AGN clus- 
tering trends observed, assuming the validity of the a x 
connection, i.e., that over a wide range of X-ray and op- 
tical luminosity, both luminosities are closely connected 
independent of redshift: 

1) The dominant accretion processes that are trigger- 
ing AGN activity could be different at different redshifts 
and/or halo masses. For example, one possible scenario is 
that at high redshift AGNs may be triggered by interac- 
tions or mergers between galaxies in group environments, 
while at low redshifts internal processes might be respon- 
sible. Consequently, due to different underlying physics 
triggering AGN activity, the observed AGN clustering 
properties at low and high redshifts could differ. 

2) Alternatively, the underlying physics may not evolve 
with time, but the kinds of AGNs being compared at dif- 
ferent redshifts and luminosities are fundamentally dif- 
ferent. It may be that high luminosity broad-line AGNs 
generally reside in lower mass dark matter halos (at least 
out to z ~ 1), and within this population the brighter 
AGNs are more clustered. On the other hand, the lower 
luminosity AGNs probed in the deep z ~ 1 surveys, 
which contain a mix of broad-line and narrow-line AGNs, 
may reside in more massive halos. In this scenario, dif- 
ferent kinds of AGNs can have different triggering mech- 
anisms. The change in the luminosity dependence of 
the clustering signal seen with cosmic time is simply ex- 
plained by observing different kinds of AGNs in the low 
and high redshift universe with current data sets. 

The present clustering measurements do not allow us to 
accept or reject these hypotheses; smaller uncertainties 
on the observed clustering and detailed model predictions 
are needed. However, it is clear from the above results 
and discussion that there is not a simple relation between 
luminosity and clustering amplitude for AGNs. 

4.4. Comparison to Galaxy Clustering 

Large-scale structure studies have established the de- 
pendence of galaxy clustering on morphological type, 
luminosity, color, spectral type, and stellar mass (e.g., 
iNorberg etaL1l2002t IMa dgwick et al. 200 j IZehavi et al.1 
I2005U iMeneux et al.l 120061 120091 : ICoil et al.l l2008t ). 
These quantities are strongly correlated with each other 
and with clustering amplitude. For a comparison with 
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our AGN clustering measurements, we focus on the color- 
dependence of galaxy clustering. The last section of Ta- 
ble 0] lists a representative sample of the color-dependent 
clustering measurements of galaxies to z = 1, where we 
hav e listed results for gala x ies with L ~ L * , whe re we 
use iBlanton et al.1 (l2003bl) : iNorberg et all (|2002t ). and 
iWillmer et all (120061)" to determine L* for each sample. 
The iMeneux et all <| 20061 ) results in Table H correspond 
to galaxies 1 mag below M* but contain more objects and 
therefore yield more reliable clustering measurements. 

Blue galaxies show a weaker clustering signal than 
red galaxies over the wide redshift range studied here 
(0.07 < z < 1.0). At low redshifts, results from larger 
surveys agree with each other and indicate that red 
galaxies are hosted by dark matter halos of slightly higher 
mass than Mdmh ~ 10 13 /i" 1 M Q , while blue galaxies 
reside in halos of M D mh ~ 10 12 h~ x M Q . At z > 0.5 
different studies find that red galaxies are more strongly 
clustered than blue galaxies, but the published cluster- 
ing scale lengths are discrepant. In terms of inferred dark 
mat ter halo masses, o ur measurements ar e in agreement 
withlCoil et al.l (|2007l ). while the results bv lMeneux et al.l 
(2006) give dark matter halo masses which are too low 
when evolve into low- 2; L* galaxies. The errors on the 
measurements of optical AGN clustering are too large to 
constrain the host galaxy type using the clustering mea- 
surements; the results are consistent with optical AGNs 
being hosted by either blue or red galaxies. However, 
iCoil et al.l (|2007l ) find that at the 2a level QSOs are clus- 
tered like blue galaxies, not red, at z = 1. 

Our clustering measurements strongly exclude the pos- 
sibility that RASS-AGNs reside in the very massive dark 
matter halos that host LRGs (Figure ITT]) . Interestingly, 
our sample of ~1000 low Lx RASS-AGNs clusters sim- 
ilarly to blue galaxies, while the high Lx RASS-AGNs 
(n = 562) have the same clustering amplitude as red 
galaxies at similar redshifts. This suggests that RASS- 
AGNs are hosted in a mix of higher fraction of blue galax- 
ies and lower fraction of red galaxies. 

5. CONCLUSIONS 

Using cross-correlation measurements between RASS- 
AGNs and SDSS LRGs, we measure the projected cross- 
correlation w p (r p ), and from this we derive the real-space 
auto-correlation function £(r) of X-ray selected AGNs 
with broad optical emission lines at (z)= 0.25 with an 
unprecedented precisio n. Using the RASS/SD SS cross- 
identification sample bv I Anderson et all ()2007l ) results in 
the largest sample of X-ray selected AGNs (n — 1552) 
used for clustering studies, covering an area of 5468 deg 2 . 
We detect a clustering signal at the ^llcr level and find 
a correlation length of r — 4.28^0.54 h^ 1 Mpc and 7 = 

1.67±o:il, fitting on scales of r p = 0.3 - 15 h^ 1 Mpc. 

We investigate the luminosity dependence of cluster- 
ing using low and high X-ray luminosity subsamples 
defined using the common AGN/QSO dividing line of 
L £5-iokeV = 10 44 ergs -i_ We detect an x _ ray luminos- 
ity dependence of the clustering signal at the ~2.5<r level, 
with the brighter sample being more clustered. This is 
contrasted with clustering measurements of low luminos- 
ity X-ray AGNs and optically selected QSOs at redshift 
z > 0.5, which suggest that low X-ray luminosity AGNs 
are more strongly clustered and reside in more massive 
dark matter halos than their high X-ray luminosity coun- 



terparts. Possible explanations are that (1) different 
mechanisms trigger AGN activities at different redshifts 
and/or halo masses or (2) at all redshift, low luminos- 
ity AGNs and the brightest QSOs reside in red galax- 
ies, while most intermediate luminosity AGNs/ QSOs are 
hosted in blue galaxies. 

Our low Lx RASS-AGN sample exhibits a similar clus- 
tering amplitude as blue, star-forming galaxies at similar 
redshifts, while the high Lx RASS-AGN sample clusters 
like red galaxies. The total RASS-AGN sample is likely 
dominated by blue host galaxies but includes a fair frac- 
tion of red host galaxies. 

We show that the auto-correlation function derived 
from cross-correlation measurements of ^1500 AGNs and 
~50000 LRGs in an area of ^6000 deg 2 constrains the 
auto-correlation function at the 10% level. Although the 
RASS is the most sensitive X-ray survey ever performed, 
only the most luminous and X-ray unabsorbed (soft) 
AGNs w ere detected. The u pcoming all-sky eROSITA 
mission (jPredehl et al.1 120071 ) with its XMM-Newton- 
like soft and hard energy high sensitivity should de- 
tect ^200,000 AGNs with much lower X-ray luminosi- 
ties. The data generated by the mission will allow one 
to compute clustering measurements with a significantly 
higher accuracy. 

To draw meaningful conclusions from AGN clustering 
measurements using AGN samples selected at different 
wavelengths and to address the evolution and luminosity 
dependence of AGN clustering, lower uncertainties on the 
detected clustering signal are required. We describe and 
successfully show in this paper that the cross-correlation 
function can be used to not only measure the relative 
clustering of AGN/galaxy samples but also to determine 
the auto-correlation function with relatively low uncer- 
tainties. 
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